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We perform a comprehensive SC/ (3)-flavor analysis of charmed mesons decaying to two pseu- 
doscalar SU(3)-octet mesons. Taking into account S'(7(3)-breaking effects induced by the splitting 
of the quark masses, m s 7^ rriu,d, we find that existing data can be described by S'l7(3)-breaking 
of the order 30%. The requisite penguin enhancement to accommodate all data on CP violation 
tends to be even larger than the one extracted from AaQp(K + K~ , ty + tv~) alone, strengthening 
explanations beyond the standard model. Despite the large number of matrix elements, correlations 
between CP asymmetries allow potentially to differentiate between different scenarios for the un- 
derlying dynamics, as well as between the standard model and various extensions characterized by 
SU(3) symmetry and its subgroups. We investigate how improved measurements of the direct CP 
asymmetries in singly-Cabibbo-suppressed decays can further substantiate the interpretation of the 
data. We show that particularly informative are the asymmetries in D — > tt + tt~ versus D — > K + K~ , 
D s -> K s n+ versus D+ K S K+, D+ -» tt+tt , D -> ttV, and D -> K S K S . 



I. INTRODUCTION 

The recent measurements of direct CP violation in non- 
leptonic charm decays [IH2] have been among the most 
exciting results in flavor physics in recent years. The 
combined significance of CP violation in charm decays is 
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for a decay d of a C = +1 meson, where A(d) de- 
notes the weak decay amplitude of the flavor eigenstates. 
This measurement, together with the plethora of avail- 
able charm data, makes an 5C/(3)-flavor symmetry anal- 
ysis worthwhile, aiming at an understanding within or 
beyond the standard model (SM) . The basics of such an 
analysis were laid out some time ago [5HH]; however, the 
present situation allows for a much more complete anal- 
ysis than was possible before. 

The SU(3) symmetry is known to be broken rather 
severely in charm decays, the most striking example given 
by T(D° K+K-)/T(D -> tt+tt-) - 2.8. However, as 
first pointed out in |10j . this does not necessarily imply 
5(7(3) breaking on the amplitude level beyond the ex- 
pected order of ~ 30% (see also [IT] for a similar analysis 
in the diagrammatic approach). Indeed, for a subset of 
decays it has been shown again recently that "nominal" 
5(7(3) breaking is sufficient to explain this ratio [T2Tll4j . 

Here we address the following questions: 



i How large is the requisite SU(3) breaking in charm 
decays? 

ii How large is the requisite penguin, i.e., triplet ma- 
trix element enhancement to explain the observed 
CP violation? 

Hi Can we distinguish between new physics (NP) con- 
tributing to operators in different representations 
of 5C/(3)? Which measurements would be particu- 
larly useful? 

In our analysis we take into account SU (3)-breaking 
induced by the splitting in the quark masses m s 7^ m u ,4 
to first order; we compare our findings to the most com- 
plete, relevant data set of decays into two octet pseu- 
doscalars and employ no further dynamical assumptions, 
the combination of which is where we go beyond previ- 
ous works on SM CP violation, e.g., recently [T2TU9| . be- 
sides a detailed assessment of the 5?7(3)-flavor anatomy 
of non-leptonic charm decays. 

The plan of the paper is as follows: The 5C/(3)- 
structure of hadronic two-body charm decays is given 
in Section |Hj In Section HI we present fits to the data 
assuming the SM. In particular, all CP violation is in- 
duced by the Cabibbo-Kobayashi-Maskawa (CKM) mix- 
ing matrix. This analysis hence holds more generally in 
all models with this minimally flavor violating (MFV)- 
feature. In Section IIVI we allow for CP violation from 
NP characterized by different SU(3) representations. We 
identify patterns in observables that can guide towards an 
identification of the underlying flavor dynamics. We con- 
clude in Section|Vj In several appendices we give Clebsch- 
Gordan tables and subsidiary information. 
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II. SU(3)-DECOMPOSITION 

We present the 5?7 (3)-decomposition of various two- 
body charm decay amplitudes. The requisite Clebsch- 
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Gordan coefficients are obtained using the tables and the 
program from Refs. [20l422j . 



A. SU(3)-limit 

SU (3)-flavor symmetry allows to express the ampli- 
tudes of the various decays d in terms of reduced matrix 
elements A fc as 



where 



A Q (d) = S^c d;lfc Af , 
A Q (d) = V* s V ud Cd-, lk A k l , 

t .k 

A (d) = V* d V us J2 c d;ikA^ , 
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We employ the commonly used classification for non- 
leptonic two-body charm decays according to the CKM 
hierarchy of their decay amplitudes: Cabibbo-favored 
(CF) at order one, singly-Cabibbo-suppressed (SCS) at 
order A and doubly-Cabibbo-suppressed (DCS) at order 
A 2 in the Wolfenstein expansion, where A ~ 0.2. 

In Eq. ([3]) i, k label the representation of the final state 
and the ucqq' interaction Hamiltonian, respectively. The 
relevant tensor product of the latter is written as 



3 2 ©6©15. 
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The initial D-meson (D = (D°,D+,D S )) is an SU(3) 
anti-triplet. The two pseudoscalar octets in the final 
state can be decomposed as (8 © 8) s = 1 © 8 © 27, 
where we symmetrized to account for Bose statistics and 
removed the dependence on the order of the final state 
mesons. The resulting coefficients c d -ij are given in Ta- 
ble|lj where we introduced A = (V* s V us +V* d V ud )/(2T,) ~ 
A 4 ~ 10~ 3 , which characterizes the CKM suppression of 
direct CP violation in D decays. Note furthermore that 
we combined the coefficients of the 3i and 32, as they 
have identical quantum numbers and therefore enter each 
amplitude with identical relative weight. Our findings are 
in agreement with Ref . [5] ; we disagree with a recent cal- 
culation [T2] in the sign of the D + — > K + ir° amplitude 1 . 
The inclusion of decays into pseudoscalar singlets is left 
for future work |2"3"1. 
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TABLE I: The coefficients Cd-ij of the decomposition into re- 
duced matrix elements in the SU(3)-limit given in Eq. |3}. 



B. Breaking SU(3) 

Including 5[/(3)-breaking through m s ^ m u ^ d , i.e., 
leaving isospin symmetry intact, leads to SU(3) break- 
ing by a single representation (8) , see also [TUl [321 [24] . 
The effective Hamiltonian contains at lowest order in the 
breaking the following decompositions: 



15 © 8 = 42 © 24 © 15i © 15 2 
6©8 = 24©15©6©3, 
3©8 = 15©6©3. 



15' © 6 © 3 , 



(6) 



We consider the CKM-leading terms in the flavor- 
breaking only 2 , hence the 3 © 8 does not contribute. A 
notable exception is the decay D + —> tt°tt + , where the 
omission of A would induce a spurious CP asymmetry. 
However, in the absence of new isospin violating interac- 
tions, the latter vanishes, see, e.g., 25, 26J. Furthermore, 
there are no non-zero matrix elements with the 15'. The 



The authors of Ref. 1 12| informed us that they agree with our 

expressions. The apparent sign differences are due to typos in 2 One should revisit this assumption once CP data become more 
their coefficient tables. precise. 
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decay amplitudes can then be written as 

A(d)=A (d)+A x (d), (7) 

where the Ax denote the flavor-breaking contributions. 
We express them in terms of reduced matrix elements 
B d 



A x (d) = Zj2 c dM B i 



A x (d) = V: d V us ^2c d;ij Bi. 
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analogous to Eq. 

Our findings for the coefficients Cd-ij are given in Ta- 
ble IV We confirm the results given in Ref. [T2] except 



for the sign of the E-terms in the SU (3)-breaking part 
of the D — > 7r°7T amplitude and the overall sign of the 
D + -> K + ir° amplitude 3 . 

For A — > 0, the flavor structure of Eq. ^ leads to 
15 SU (3)-breaking matrix elements in addition to the 
three from the SU(3)-limit. However, the coefficient ma- 
trix of Clebsch-Gordan coefficients does not have full 
rank, implying that not all matrix elements arc phys- 
ical. Since in fact the coefficient matrix has rank 11, 
we go on and reduce the number of matrix elements by 
7. Specifically the two 3 representations have again the 
same quantum numbers, allowing us to combine the cor- 
responding coefficients. Using Gaussian elimination we 
can further remove fif 2 , B^ 3 , B^ 3 , B^ 2 , and Bfj. We 
do not mix leading SU(3) elements in the process, and 
keep their normalizations when subleading matrix ele- 
ments are absorbed. The A-suppressed S'J7(3)-breaking 
matrix elements that appear in the course of the redef- 
initions are again neglected as they are of higher or- 
der in the power counting employed in this work. Due 
to the reparametrizations, the sub- and superscripts on 
the B\ cease to correspond to the SU (3)-representations 
of the interaction and the final states. The appear- 
ing combinations have been normalized according to 
|ci| 2 -B phys = 'Yli i CiB^ UtyS \ in order not to change 
the relative normalization of the matrix elements. The 
resulting coefficients Cd-ij of the physical decomposition 
are given in Table [V] As anticipated, we end up with in 
total 13 independent matrix elements. In the following 
analyses we use Eqs. (|8|-(10) with the coefficients from 
Table M 



C. SU(3)-breaking resistent sum rules 

For A = the SCS 8x11 submatrix of Clebsch-Gordan 
coefficients including 577(3) breaking in the physical 



parametrization has rank 7. Therefore, we know a priori 
that there is only one linear sum rule among the SCS am- 
plitudes that remains valid after SU(3) breaking. This 
sum rule is the well-known isospin relation [27] 4 



1 

71 



A(D° -> ir+ir-)+A{D° -> 7r°7r°) = A{D+ -> tt+tt ) . 

(11) 

One may ask whether there are also approximate sum 
rules which are broken by a single matrix element only. 
By calculating the rank of the corresponding matrices of 
Clebsch-Gordan coefficients, we find that among the SCS 
decays there is exactly one such sum rule. It is given by 
the (quasi-)triangle relation 

A(D+ -> tt+tt ) - -^=A(D S -> K°tt+)- 
V 2 



A(D s ^K+tt q ) = Z x I-B< 7 



(12) 



which generalizes Ref. 27J, where SU(3) breaking by 
triplet matrix elements only was considered. Setting 
B27 1 = l ne sides \Ai\ of the triangle have to obey 
|A max | — 1-^2 1 — I -A3 1 < 0, i.e., the longest of the three sides 
has to be smaller than the sum of the other two. This re- 
lation, if broken, would prove the necessity for B^ 1 7^ 0; 
it is, however, fulfilled by the data, given in Table [TTJ. In 
the sys tem before reparametrizations the right-hand side 



of Eq. (12 1 involves two matrix elements, B^ 1 and £? 4 f , 



and there is no sum rule with just one breaking matrix 
element. 

In case the matrix element B^ 1 can be neglected and 
assuming MFV/SM, Eq. (12) allows to extract the size 



of the penguin contributions by measuring the involved 
branching ratios and CP-asymmetries by the usual tri- 
angle construction. It involves a common base |.4(D + — > 
7r+7r | = \A{D~ -> 7T-7r°)| for the triangle and its CP- 
conjugate one; the requisite weak phase 7 can be taken 
from global CKM fits. A finite B^ 1 induces a correc- 
tion to this procedure of the order of the S'C/(3)-breaking, 
which however, is not so small. 

Note that once CF and DCS modes are considered as 
well, further sum rules arise, see 1141. 



III. SM/MFV FITS 

We confront the SU (3)-analysis from the previous sec- 
tion to data. The relevant measurements are compiled in 
Table [11} The fits are carried out using the augmented la- 
grangian [2H1 [25] and Sbplx/Subplex algorithms (SU EI] 
that are implemented in the "NLopt" code |30j . 



3 See footnote [TJ 



The apparent discrepancy between Eq. |TTJ) and previous 
works [9] 127] is due to the different conventions used, specif- 
ically A (D° -> 7r°7r°) here 

A (D+ -> TT + TT ) here = -A (D+ -> TT+TT ) 



V2A (D° -> tt°tt°) and 

/ previous 
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Our goal is to see whether SU(3) gives a reason- 
able expansion for the full set of two-body decays of 
charm to pseudoscalar octet mesons, as shown possible 
for D° -> P+P- [iniCEMI]- Our framework is MFV, 
which includes the SM. In these models CP violation is 
suppressed by A, the relative weak phase is order one, 

&r g (v; b v ub /v c * d v ud ) = - 7 . 

To obtain the sizable CP violation observed in SCS de- 
cays, the triplet matrix elements, A\ or A|, need to be 
sufficiently large, as they are the only ones not severely 
restricted by the branching ratio data. In terms of the 
low energy effective theory, this concerns penguin con- 
tractions of tree operators, the chromomagnetic dipolc 
operator ua^G^c, where G^ v denotes the gluon field 
strength tensor, and the QCD penguins uj^cj^q 97^9; 
all of which are purely SU(3) triplets; we neglect the con- 
tributions from electroweak penguin operators. On the 
other hand, matrix elements involving the 6 and 15 rep- 
resentations receive contributions from the CKM-leading 
tree operators only. 

Given the lack of a dynamical theory for hadronic 
charm decays, we have to resort to order-of-magnitude 
arguments when judging the fit results. The penguins 
(triplet matrix elements) are generically expected to be 
suppressed by a factor of ~ a s /Tr ~ 0.1 compared to 
their tree counterparts. While this estimate cannot be 
expected to hold literally, in the past an upper bound 
of one for this ratio was widely considered conservative, 
leading to rather strong upper limits for SM CP viola- 
tion. Possible enhancements have been discussed from 
an early stage on [221 133], and recently, e.g., in [i"4l IT9] . 
Note that the widely used analogy to the AI — 1/2 rule 
in kaon decays seems questionable, as the effect is ex- 
pected to scale as m s /m c |32j . However, enhancements, 
e.g., from rescattering cannot be excluded, but the nec- 
essary size for reaching the present central value seems a 
stretch, even in these analyses. Below, we introduce mea- 
sures to quantify penguin enhancement for our analysis, 
and discuss the results obtained with present data. 



A. Observables vs. degrees of freedom 

The branching ratio of a decay d into two pseudoscalars 
Pi, 2 in terms of its amplitude A(d) is given as 

B(D^ P 1 P 2 ) = T D V(d)\A(d)\ 2 , (13) 

with the phase space and normalization factor 

a/ (m 2 D - (mi - m 2 ) 2 ){m 2 D - {mi + m 2 ) 2 ) 



We further employ differences and sums of direct CP 
asymmetries of SCS decays to final states /i 2, defined as 



V(d) 



16 7rm|j 



(14) 

We take the lifetimes r^. and masses m, of the involved 
mesons from [3~i] . 

For decays involving D° or K°, the direct CP asymme- 
try defined in Eq. ^ is not measured directly in the ex- 
periment. The corresponding indirect contributions are 
subtracted before fitting, see Appendix [X] for details. 
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respectively. Considering Aa^'p(/i,/ 2 ) instead of the 
individual asymmetries is experimentally advantageous; 
furthermore, to a very good accuracy, indirect contribu- 
tions cancel in the difference. 

We consider now the parameter budget of the SU(3) 
ansatz including breaking effects at leading order versus 
the available experimental data. The 17 decay modes 
correspond in principle to 26 observables, i.e., 17 branch- 
ing ratios, 8 (direct) CP asymmetries of SCS decays, 
and the strong phase difference between D° — > tt + K~ 
and D° — > ir~ K + . Since, however, the branching ratio 
B(D S — > KlK + ) is not measured yet, we are left with 25 
observables. 

On the parameter side within MFV/SM, there are 
three S'J7(3)-limit matrix elements that come with a fac- 
tor of £ and two S , J7(3)-limit matrix elements that come 
with a factor of A only. Of the five resulting matrix 
elements one can be chosen real as we are only sensi- 
tive to differences of strong phases. In the SU (3)-limit 
there are hence nine real parameters. Taking into account 
5?7(3)-breaking, we end up with 13 independent matrix 
elements, see Section II B i.e., 25 real parameters. 



B. The fate of unbroken SU(3) 

The need for flavor-breaking in SCS decays is most 
obvious in the large difference in the rates of the K + K~ 
and 7r + 7r _ modes, and in the enhanced branching ratio 
for D° — > KgKg, whose contribution oc S vanishes in 
the SU{3) limit. Fitting in this limit CF and DCS modes 
only also returns a very large % 2 , indicating sizable flavor 
breaking in these modes as well. This is seen, e.g., in 
B(D° -> K+tt-)/[\ 4 B(D° -> K-TT+)] - 1.5. 

Fits in unbroken SU(3) using CF modes alone are not 
possible without additional assumptions, as the respec- 
tive decay amplitudes have too much of a linear depen- 
dence, see also the observations made in [35] - One may 
entertain the possibility of a CF-only fit by including 
5J7(3)-singlets in the final states. While this leads to ad- 
ditional observables and constraints, it also leads to ad- 
ditional matrix elements. The fit presented in Ref. [36] 
yields x 2 — 1-79 for 1 degree of freedom (dof) at the price 
of an additional dynamical assumption, which effectively 
relates matrix elements involving singlets and octets. 



C. Fitting flavor-breaking 

Including linear S'J7(3)-breaking, we generically obtain 
good fits for a multitude of configurations. We find that a 
reasonable x 2 /dof requires at least two S'J7(3)-breaking 
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matrix elements present. However, a fit with the as- 
sumption of triplet enhancement, as proposed in [27|l32j. 
and recently investigated in [T2], does not yield accept- 
able results {x 2 /dof — 8.6), as no SU(3) breaking en- 
ters the CF/DCS sector, see Table [V] Therefore, at least 
one matrix element other than Bf 8 is needed to describe 

the data; using, for instance, Bf and B^ 2 we obtain 
X 2 /dof = 1.3. Note that a fit with Bf to the SCS modes 
alone does work, x 2 /dof = 1.0. 

To evaluate the convergence of the S't/(3)-expansion, 
we quantify the size of the flavor-breaking with the fol- 
lowing measures: The size of the SU (3)-breaking matrix 
elements, defined as 



"A" 



maxy|^| 



(17) 



max(|^|,|4|,|Af|) ! 
and the size of the SU (3)-breaking amplitude, written as 



S' x = max d 



A x (d) 



A(d) 



(18) 



To avoid a bias in the latter definition, we exclude the 
decay D° — > K$Ks for which the amplitude in the SU (3) 
limit is a A. We use both measures, as S x ignores 
the possibility of a suppression of ST/(3)-breaking from 
Clebsch-Gordan coefficients, while 6' x ignores the possi- 
bility of large cancellations. 

We find that the data can be described by an SU (3)- 
expansion with 5 X ] < 30%, see Fig. 1 where we show 
68% (dark red) and 95% (light orange) confidence level 
(C.L.) contours relative to the best fit point, see Sec- 
tion HID While solutions with larger SU(3) breaking 




FIG. 1: The 68% (dark red) and 95% (light orange) C.L. 
contours in the Sx-8' x plane with respect to the best fit point. 



are not excluded by the data, we take this result as con- 
firmation that the expansion works as good as could be 
expected. We therefore assume its validity in the follow- 
ing, and exclude solutions with cancellations by imposing 
6% < 50% in most fits. 



Having established that with the current data the 
5J7(3)-expansion can be applied, we study the anatomy 

of the fit solutions. With the upper bound on 5 X , 
we obtain a fairly good fit to the full data set with 
X 2 /dof = 1.6 for the configuration mentioned above, 
which consists of only the two matrix elements Bf and 
B^ 2 . We stress that the minimal number of Bf is two. 
A nicer fit, however, is obtained if at least three flavor- 
breaking matrix elements are present. Such examples 
are the configurations with Bf, B 1 ^ 2 , B^ 1 or Bf, B^ 1 , 
-B27 2 1 both of which give \ 2 /dof — 1.0. In the following 
we do not look at specific configurations with a minimal 
number of flavor breaking matrix elements but rather fit 
the full system with all Bf . 



D. Penguin enhancement 

Next we turn to analyze the penguin enhancement. 
One possible definition is the following ratio: 



<\a%\M\) 



max(\A%\,\A*\,\A}*\) 



(19) 



In order to inspect also here the possibility of cancella- 
tions, we further define the penguin-to-tree fraction of a 
specific decay d as 



S' 3 (d) = 



Cd:27 15^27 + C d;& 6 A t + c d;S 15^ 

and its maximum 

5' 3 = maxd J 3 (d) . 



(20) 



(21) 



Here, for the same reason as in the case of S' x , the decay 
D° —> K$Ks is not taken into account. 

At the best fit point we obtain \ 2 = 1-0 for 25 real fit 
parameters and an equal number of observables, obtained 
without a constraint on S x \ The residual \ 2 is caused 
by Acp{D + — > 7r + 7T°), which is measured nonzero at la, 
but impossible to accommodate in the SM as it vanishes 
therein. We observe that 8<p is driven to huge values of 
0(100) in the fit. The reason for this strong enhance- 
ment lies not so much in the measurement for Aa^p, 
but is due to the CP asymmetries A C p{D° — > KsK s ), 
A C p{D s -> K s ir + ), and A C p{D s -> A'+rr ), which 
presently have very large central values, see Table [H] 
Their uncertainties are very large as well, rendering the 
effect insignificant for each single measurement. Together 

however, while still allowing for solutions with sj[ ~ 5 
at 95% C.L., they shift the 68% C.L. region to very large 
values. Numerically, we obtain x 2 = 1-9 for ^ < 30, and 
X 2 = 3.6 for 8z < 10, S3 in both cases saturating the 
bound. The latter value is slightly increased to % 2 = 4.5 
when additionally imposing 5 X < 50%, indicating a very 
small correlation between the two measures. These ob- 
servations are illustrated in Figs. [2] and [3] In the latter 



G 




Observable Measurement References 



SCS CP asymmetries 



Aa^ T p (K+K- ,7r+7r-) 


-0.00678 ± 0.00147 


[Ml 1371138] 


T,a^ v p (K+K-,ir+TT-) 


+0.0014 ± 0.0039 


tpMl 1371139] 


A CP {D° 


-» A"s^s) 


-0.23 ±0.19 


m 


A CP (D° 7r°7r°) 


+0.001 ± 0.048 


HQ] 


A CP (D+ 


-> 7T°7r+) 


+0.029 + 0.029 




A CP (D+ 


-> A S A+) 


-0.0011 ± 0.0025 




A C p(D s 


— > A' S 7T+) 


+0.031 ± 0.015 


t gum] mi 


A C p(D s 


-> A+7T ) 


+0.266 + 0.228 




Indirect CP Violation 


a 


nd 
CP 


(-0.027 + 0.163) ■ 10~ 2 


® 


5 L = 2Rc(e)/(l + \e\ 2 ) 


(3.32 + 0.06) ■ 10" 3 


El 




K+w- 


strong phase difference 




5 




21.4° ± 10.4° 


t0 


SCS branching ratios 


B(D° - 


K+K-) 


(3.96 + 0.08) ■ 10~ 3 


E3 


B(D° - 


7T^~ TV ) 


(1.401 ± 0.027) ■ 10" 3 


EH 


B{D° - 


> ^s^s) 


(0.17 + 0.04) ■ 10- 3 


El 


B(D° - 


-> 7r°7r°) 


(0.80 + 0.05) ■ 10~ 3 


34 


B(D+ - 


-> 7r°7T+) 


(1.19 + 0.06) ■ 10~ 3 


El 


B(D+ - 


> A S A+) 


(2.83 + 0.16) ■ 10~ 3 


ED 


B(D S - 


* A 5 7T+) 


(1.21 ± 0.08) ■ 10" 3 


EH 


B(D S - 


* A+7T ) 


(0.62 + 0.21) ■ 10~ 3 


EH 


CF* branching ratios 


B(D° - 


> K~TT+) 


(3.88 + 0.05) ■ 10~ 2 


EH 


B(D° - 


-> A S 7T°) 


(1.19 + 0.04) ■ 10~ 2 


EH 


B(D° - 




(1.00 + 0.07) ■ 10~ 2 


El 


B(D+ -> K s n+) 


(1.47 + 0.07) ■ 10~ 2 


El 


B(D+ -* K l tt+) 


(1.46 + 0.05) ■ 10~ 2 


El 


B(D a - 


► K S K+) 


(1.45 + 0.05) ■ 10~ 2 


tElHH] 


DCS branching ratios 


B(D° - 


> K+TT-) 


(1.47 + 0.07) ■ 10~ 4 


El 


B{D+ -> K+n°) 


(1.83 + 0.26) ■ 1GT 4 


El 



TABLE II: The observables and the data for indirect CP vi- 
olation used in this work, see Appendix [X] for removal of ef- 
fects from charm and kaon mixing. 'The measurement quoted 
corresponds to our average. Systematic and statistical uncer- 
tainties are added in quadrature. '"Our symmetrization of 
uncertainties. *Modes into Ks,l assigned to CF decays. 



we show the influence of the largish measured CP asym- 
metries explicitly by excluding them from the fit. As a 

result, values of 5^ ~ 3 become allowed at 68% C.L., 
consistent with Refs. |13l [14] [47] , where these asymme- 
tries have not been taken into account either. 

Closer inspection exhibits common features among the 
presented fits: The penguin matrix elements tend to 
be largely enhanced; there are large hierarchies between 
5' 3 (<i), depending on whether both triplet matrix elements 
A\ and A\ are present, as is the case in the decays 
D° -+ K+K~, D° -+ 7T+7T- and D° -+ 7r°7r°, or not. In 
the latter modes indeed cancellations take place, yield- 



15 



, cn 10 



5 





10 20 30 40 50 

5> 3 



FIG. 2: The 68% (dark red) and 95% (light orange) C.L. 
contours in the 5' S -S3 plane with respect to the best fit point. 




S' 3 

FIG. 3: The same as in Fig. [5] without the data on A C p(D° 
K S K S ), A C p(D a -+ A S 7T+) and A C p(D s -> A'+TT ). 

ing a smaller penguin amplitude with S' 3 (d) -C S3, and 
typically 6' 3 (d) ~ 0(5). 

The maximum value of S 3 (d), on the other hand, usu- 
ally comes from the decays where only Ag is present, i.e., 
D+ -+ K S K+, D s -+ K s tt+ or D a -4 K+n°. Since fur- 
thermore typically max(|A$|, \Al 5 \) = \A\ 5 \, and 
the ratio of the Clebsch-Gordan coefficients of A| and 

Ag 5 in the latter modes equals ~ 3, we understand 
why generically S' 3 ^ 3 ■ S3. 

Taking into account the full dataset, we therefore find 
indications of even stronger enhanced penguin ampli- 
tudes than previous analyses relying on Aa^p only. On 
the other hand, given the present uncertainties, the large 
central values responsible for this additional enhance- 
ment might be assigned to experimental fluctuations. As 
the enhancement indicated by Aa^'p is already very dif- 
ficult to explain within the SM, any additional enhance- 
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ment challenges this option further, making more precise 
measurements of the corresponding CP asymmetries ex- 
tremely important. 



IV. PATTERNS OF NEW PHYSICS 

In the absence of a dynamical theory or further in- 
put we still cannot make a definite statement on whether 



the enhanced penguins observed in Section III stem from 
physics within or beyond the SM. Note that knowing the 
weak phase 7 precisely in MFV/SM is of no help as there 
is an approximate parametrization invariance: because of 
the smallness of |A|, the rate measurements are not sensi- 
tive to the terms 0(ReA, |A| 2 ). The CP asymmetries are 
proportional to ImA, however, they always involve un- 
known matrix elements not constrained by the rates. As 
a result, a shift in ImA (e.g., a different weak phase) can 
in the fit always be absorbed by a shift in the magnitude 
of the respective matrix elements, as long as extreme 
values are avoided. Therefore, we set in the following 
A = Asm, even when the NP model does not require a 
specific value. The potential enhancement factors from 
the Wilson coefficients are then identified by enhanced 
NP matrix elements, where we assume the absence of 
very large enhancement factors from QCD. Furthermore, 
we assume the validity of the SU (3) expansion, using 

again the upper limit 5% < 50%. 

To investigate the interplay of NP and SU(3) breaking, 
we study the following S'J7(3)-limit relations involving 
CP asymmetries: 



T(D° 



K+K- 



a d j P (D° 



T(D° 
T(D+ 



"'CP 
dir 

a CP 



K+K- 



a d S r P {D° -> K°K°) = 0, 



(D+ K°K+) 



"CP 



(D+ -> 7T+7T ) = 0. 



(22) 

(23) 

(24) 
(25) 



The first three relations follow in the J7-spin limit, while 
the last one is an isospin relation. All relations can be 
obtained by inspecting, for instance, Table [I] Eqs. ( 22 ) 
and ( 23 1 can be rewritten as 



Ea^piK+R-^+TT-) = 0, 
Y,a%p(K K + ,K°Tr + ) =0, 



(26) 
(27) 



with corrections of the order O(ReAlmA), which are 
completely negligible compared to SU (3)-breaking cor- 
rections. In addition, the S'C/(3)-breaking representa- 
tions 6 and 24 do not yield corrections to these relations, 
either. 

The ST/(3)-expansion allows to quantify corrections to 
the above relations. Further deviations then indicate the 
presence of U- and isospin changing interactions beyond 
the SM, models of which are discussed in the next section. 



A. New physics scenarios 

We recall the 5?7(3)-limit Hamiltonian "Hf£ s for SCS 
decays within MFV/SM: 



1/SCS 
n SM 



nl GS = E I - 



-T/SCS 
n A,SM 



-T/SCS 



73 



15 



3/2 



15i/ 2 - 6 X /2 



(28) 



and H^sm involving the representations 3 and 15, 
which, however, do not contribute significantly due to 
our assumption of "well-behaved" hadronic matrix ele- 
ments, and are therefore neglected in the following. Here 
and in the following the subscript to the representation 
labels the shift in total isospin AI. 

We discuss scenarios Ujfp = Uf^ps + ^| CS > which 
can be classified according to the contributing SU(3)- 
representations. We use A = AE = (V* s V us + V* d V ud )/2. 

Very similar to the SM/MFV are models based on new 
physics in the 3 alone ("triplet model"), 



u 



SCS 



A 



V 



,NP 



(29) 



At the quark level, the corresponding operators are the 
chromomagnetic dipole operator and the QCD penguin 
operators. They have been discussed in supersymmetric 
and extradimensional models [5T] in the context 
of CP violation in charm. 

We further consider a 3 + 15 interaction 



W|+15 = A 



15 



NP 

3/2 



1 

T2 



15 



NP 
1/2 



,NP 
'1/2 



(30) 



which arises from an operator with flavor structure ucuu 
( "HN model" ) . Those have been investigated recently in 
a peculiar variant of a 2-Higgs doublet model (2HDM) 
that links the top sector to charm While the fla- 

vor structure is identical to the corresponding SM tree 
operator, it can have a different Dirac structure. As a 
result, the matrix elements of the 15-representation are 
independent of those appearing with a coefficient E, and 
therefore not constrained by the data on the rates. 

Finally, we allow for 3 + 15 + 6 terms ("AC/ = 1 
model"), 



n 



SCS 

3+6+15 




^NP 
} l/2 




(31) 



The corresponding operators have flavor content scus. 
This structure may arise from tree level scalar exchanges 
as in 2HDMs or a color octet, see [S3] for a list. The ap- 
pearance of a third representation carrying a weak phase 
different from the leading contributions implies signifi- 
cantly less correlation between different CP violating ob- 
servables. This makes this scenario especially difficult to 
identify in patterns. 
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All models contribute to SCS decays only. We do not 
consider potential constraints on the scenarios with NP 
from 4- Fermi operators by D° — D° mixing and e'/e [4"T] . 

As QCD preserves 5'C/(3)-flavor, the irreducible rep- 
resentations form subsets which renormalize only among 
themselves 33J. While such effects often have numerical 
relevance, they do not for our analysis, as we fit ma- 
trix elements rather then calculating them. We still ask 
whether the SU (3)-anatomy of the NP models considered 
is radiatively stable. 

A NP contribution in the 3 does not "spread out" by 
mixing onto other representations. The situation for the 
scalar operators srClUrsl is likewise simple: They mix, 
together with their scalar and tensor color-flipped part- 
ners with the same flavor content at leading order only 
among themselves. The mixing onto the dipole operators 
vanishes for m s = md = [53, 54]. The total contribu- 
tion to the amplitude therefore remains 3+15 + 6. 

The scalar urClUlUr mixes among itself, color-flipped 
partners and onto chirality-flipped QCD-penguins; 
hence, renormalization group (RG) running modifies the 
weight between the 3 and the 15 set at the NP scale by 
Clebsch-Gordan coefficients. However, an explicit cal- 
culation of the leading order running between the weak 
and the charm scale shows that this effect is < 7%, and 
can be safely neglected for the purpose of our analysis. 
Anomalous dimensions can be taken from, e.g., |55) . 

The RG stability depends in general on the Dirac struc- 
ture. We use in the following the corresponding SU (3)- 
classifications, but have in mind the model examples 
where RG effects can be neglected. 



B. Patterns 

In this section we investigate if and how the different 
NP scenarios can be distinguished from each other and 
MFV/SM. To that end, we perform fits to the full data 
set in each scenario and look for specific correlations. 

Given the fact that only one combination of CP asym- 
metries is measured significantly non-zero so far, the dif- 
ferentiation of models with present data is extremely dif- 
ficult. We therefore consider in our analysis in addition a 
future data set, see Table \UJ\ and Appendix [C| for details. 

Note that within our framework we are not able to dis- 
tinguish a NP 3 from the SM. A recent idea to resolve 
this is to measure radiative charm decays, where the sen- 
sitivity to enhanced dipole operators is enhanced w.r.t 
hadronic decays [55] . 

We start by analyzing some generic features of 
Eqs. (J22])- (|25"] ). The most clearcut relation is the isospin 
one, Eq. p5j); it holds even in the presence of SU(3) 
breaking, as long as no operator with a 15 3 / 2 represen- 
tation different from the SM one is present, see also [26] . 
Therefore, it serves as an extremely clean "smoking gun" 
signal for the HN model. 



The unique feature of Eq. ( 24 ) is the vanishing of the 
CKM-leading part of the D° — >• KgKs amplitude in 



the SU(3) limit. The size of the rate is determined by 
the SC/(3)-breaking contribution; the corresponding CP 

asymmetry is therefore enhanced at 0{A/8^). This en- 
hancement might roughly be estimated as 



"•CP 
l CP 



(D° -> K°K°) 

(D° ->• K+K-) 



IBR(D° -> K+K-) 
BR(D° -> K°K°) 



(32) 



implying a$ P (D° K°K°) ~ 1% for a$ P {D° 
K + K~) ~ AajSp/2. While this might be considered an 
upper limit for the SM, as the asymmetry in the charged 
final state is already larger than expected therein, further 
enhancements are possible within the NP scenarios. The 
sole difference between the NP models is how strongly 
a^ p {D a -> K S K S ) is correlated to other CP asymme- 
tries. In the MFV/SM and the triplet scenario it is de- 
termined by the triplet matrix elements alone, while in 
the other two NP scenarios the 15 can give a significant 
contribution as well. 



Observable 



Future data 



SCS CP asymmetries 
Aa% r P (K + K- .tt+O -0.007 + 0.0005 
Yla%{K + K-,-K + K-) -0.006 + 0.0007 
a d j P {D+ -> K S K+) -0.003 ± 0.0005 



0.0 + 0.0005 



a A j P {D 3 ->■ K+tt°) 



0.05 ± 0.0005 



K 7r strong phase difference 



21.4° ±3.8° 



TABLE III: Future data, all other values as in Table [II] The 
central values of the single CP asymmetries that correspond 
to Aa$ P and Eocf are a% P {D° -> K+K') = -0.0065 and 
a% P (D° -> 7T+7T-) = 0.0005. 



The remaining two relations Eqs. (|26j)-(27) are broken 
differently in different NP models. In MFV/SM, as well 
as in the triplet and HN scenarios, the sum of the CP 
asymmetries receives two contributions at 0(Sx )'■ one is 
proportional to Aa^p, driven by the relative rate differ- 
ence of the two modes, i.e., the contribution simply stems 
from the different normalization of the two CP asymme- 
tries. The other contribution stems from the interference 
of the SU(3) breaking part of the amplitude with the 
part oc A. The sign of the total contribution can not 
trivially be extracted from the ratio of the rates. 

In addition to these 5[/(3)-breaking contributions, in 
the AU = 1 model the relations Eqs. (26 1-(27) are broken 



at 0(1) by NP. The reason for that is that in this model 
there is no discrete U -spin symmetry of the Hamiltonian 
under exchanging all down and strange quarks anymore. 
Therefore, generically the contributions of the AU = 1 
model are expected to be larger than in the other NP 
scenarios. 
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Not unexpectedly, all scenarios fit the current data 
well, with a minimal \ 2 ~ 1- The main difference lies 
in the interpretation of the enhancements of the various 
matrix elements. In addition, as mentioned above, the 
HN model has the advantage of being able to explain 
the CP asymmetry in D° — > ir°ir + as well, leading to 
X 2 ~ 0. Excluding this measurement, all scenarios have 
a vanishing minimal \ 2 ■ 

In the HN model we get good fits with A^'^ /A^ 5 ~ 
10 which is essentially determined by the la measure- 



0.15 



ment of a c h P (D+ 



7T°7r+) 



because the A^' N is the 



only contributing NP matrix element to this mode. As 
generic size of the urClUlUr to tree enhancement we 
obtain 




-0.010 -0.005 0.000 0.005 0.010 



qSU(3) 



^3/2 .15.JVP 
10 ^27 
2 A 15 
5v^2 8 



13. 



(33) 



FIG. 4: 95% (solid) and 68% (dashed) C.L. contour lines for 
the current data with only three breaking matrix elements 
(Bf, Sf 2 , B™ 1 ) with 6$ < 50%. 



where we accounted for the Clebsch-Gordan coefficients. 
This value is in agreement with [52] . It also fits the recent 
results on the forward-backward tt production asymme- 
try from the full data set from the CDF experiment [ST] . 
The related bound from the LHC on the charge asymme- 
try Ac GUI [55] can "presumably" be evaded by the HN 
model [60]. 

We learn that with present data a clear separation be- 
tween different NP models is not possible. There are two 
paths to obtain a clearer picture. Either we gain insights 
in the dynamics of SU(3) breaking, which would, e.g., 
allow us to identify certain matrix elements as leading in 
the breaking, or we wait for more precise data to see if 
the patterns of NP discussed above become significant. 

To demonstrate how an improved understanding of 
SU(3) breaking would improve our fits, we consider one 
of the scenarios mentioned above with only three addi- 
tional matrix elements, J3f,i?g 52 , and B 2 ^ 1 . As for this 
case MFV/SM already fits the data well, so do the NP 
scenarios. We obtain x 2 / d of = 10/10 in MFV/SM and 
the triplet model, x 2 / do f = 3.2/6 for the HN model, 
and x 2 /dof = 5.4/4 for the AU — 1 model. The slightly 
worse result for the AU = 1 model is due to the fact that 
the additional degrees of freedom introduced are not nec- 
essary to fit the CP asymmetries, given that only Aa c n P 
is measured significantly different from zero. 

We find that for such an S'C/(3)-breaking scenario, 
the different NP models start to imply different pat- 
terns. This is illustrated in Fig. HI where we observe 
a correlation between the signs of a^ P (D° — > Tr°7r°) and 

a c 1 p (K + K~ , 7r + 7r~) in the triplet model. The capabil- 
ity to differentiate between the different models will im- 
prove significantly with future data. With present data, 
however, it is already possible to exclude many scenarios 
for SU(3) breaking, as the exclusion of the pure triplet 
enhancement demonstrated in Section [ill CI 

The future data scenario is designed having in mind the 
AU = 1 model being realized. The goal is to determine 
whether this model can be distinguished from the others 
by D — !• PP data, despite its many CP violating contri- 



butions. We find that all NP models remain capable of 
fitting the future data set well, despite its somewhat spe- 
cific construction, with Xmin ~ f° r * ne HN model and 
Xmin ~ 1 f° r the others. In Fig. [5J we show exemplarily 
correlations for the various models as in Fig. [4] While 
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FIG. 5: 95% (solid) and 68% (dashed) C.L. contour lines for 

the future data set with 6j) < 50%. Note that the contour for 
the AU — 1 model lies underneath the one of the HN model. 



for the shown observables the AU = 1 model cannot be 
distinguished from the HN model, in MFV/SM and the 
triplet model a clear prediction of a sizable CP asymme- 
try in D° — > ir°ir° emerges. This serves as a confirmation 
that in principle the various scenarios are distinguishable 
within our framework. However, as is obvious comparing 
Figs. [4] and [sj further dynamical input on SU(3) breaking 
would facilitate this task enormously. 
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V. CONCLUSIONS 



Acknowledgments 



We performed a comprehensive 5t/(3)-flavor analysis 
using the complete set of data on two-body decays of 
D-mesons to pseudoscalar octet mesons. The results are 
based on plain SU{3); in particular, we did not make any 
assumptions about decay topologies nor assumed factor- 
izable SU (3)-breaking. We find 

i The S'J7(3)-expansion can describe the CF, SCS 
and DCS data set with breaking of 0(30)%. 

ii 5C/(3)-breaking matrix elements involving higher 
representations cannot be neglected. 

Hi Current data imply significantly enhanced penguin 
matrix elements with respect to the non-triplet 
ones, see Figure [2[ If the measured largish CP 
asymmetries in D" — > KgKg, D s — > Ks n+ an d 
D s — > K + ir° decays are excluded from the fit the 
penguin enhancement is shifted to smaller but still 
significantly enhanced values at 0(2 — 5), see Fig- 
ure |3j favoring interpretations within NP models. 
Improved data could clarify how large the penguin 
enhancement actually needs to be. 

The model-independence of the S'J7(3)-expansion leads 
to a large number of S'L r (3)-breaking matrix elements, 
and the outcome of the fits currently does not allow 
an unambiguous interpretation regarding the underlying 
electroweak physics. In more minimal scenarios, where 
S't/(3)-breaking is limited to a smaller number of matrix 
elements, clearer patterns exist, see Figure|4] In addition, 
it is possible already with present data to differentiate 
between various SU (3)-breaking structures. 

Keeping all matrix elements, the following measure- 
ments of direct CP violation are found to be informative 
for discriminating scenarios: 



iv A breakdown at 0(1) of the relations Eqs. (26) and 
fl27| between D° -> K + K~ versus D° -> vt+tt" 
and D s — > Ksir + versus D + — > K$K + would indi- 
cate AU ^ new physics. An observation would 
support the 3 + 6 + 15 Hamiltonian. 

v An observation of a finite a d j P (D + — > 7r + 7r°) would 
signal AI = 3/2 new physics. The current la hint 
with large central value already favors the 3 + 15 
Hamiltonian, an example of which is provided by 
|52) . The operator enhancement from the current 
fit is, given the uncertainties, in the ballpark of 
what is required to explain the current anomalies 
in the top sector. 

vi We predict that a^ P (Z?° — > KsKs) is enhanced 
with respect to Aa^p, see Eq. (32). 



Future fits with improved data on charm CP violation 
[BTH64] will shed light on the origin of CP and flavor 
violation in charm. 

Note added: While this work has been completed, 
a related study appeared [55] . 
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Appendix A: Subtracting indirect CP asymmetries 

Decays with a Ks l in the final state receive an ad- 
ditional contribution to the CP asymmetry from kaon 
mixing. To isolate the direct CP asymmetry, we sub- 



tract the mixing contribution cx 8l 



2Ree 

i+M 2 



For a decay 



with a K° (K°) in the flavor final state, the mixing con- 
tribution to its CP asymmetry is given by Aq P = §l 
(A§° p = -Sl) [42 . Therefore, 

afj P (D+ -> K S K+) = A CP (D+ -> K S K+) + S L , (Al) 
a d S r P {D s -> K s tt+) = A CP (D S -> K S 7r+) - 6 L . (A2) 

As pointed out in |66j , the actual influence of kaon mixing 
depends on the experiment, due to its dependence on the 
kaon decay time. In the most recent analyses [42, 44 this 
effect is accounted for. 

To obtain aQ P (D —> 7r°7r°), we took into account the 
effect of D° — D mixing in the same way as kaon mixing 
and subtracted a™ P . 

Indirect CP violation in both kaon and charm mixing 
is discarded in Acp{D° — > KsK$), as for this mode the 
experimental uncertainties are much larger than these 
effects. 

To calculate Y,Acp{K + K~. 7r + 7r~), we average the 
data from CDF and the B-factories, taking into ac- 
count the correlations between Acp(D° —> K + K~) 
and Acp(D° — > 7r + 7r~). In order to obtain 
Y,acp(K + K~ , 7r + 7r~), we subtract the contribution from 
indirect CP violation. We find the correlation coefficient 
between Aafjp and SaJ^p to be only ~ 5%, which we can 
safely neglect in our analysis. The reason for the small 
correlation lies in the small impact of the £?-factory re- 
sults and a™ P on the average of Aa d ^ P (K + K^, 7r + 7r~). 



Appendix B: Amplitudes including SU(3) breaking 

In this appendix, the coefficient tables for the SU (3)- 
breaking parts of the amplitudes are given. In Table |IV| 
we list the coefficients for the full set of matrix elements, 
before reducing the basis to include only physical ones. 
The coefficients for the latter are given in Table |V} 
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TABLE IV: The coefficients Cd-ij of the S'?7(3)-breaking decomposition given in Eqs. |8|-(10l without reparametrizations 



Appendix C: Future data scenario 

To obtain estimates for the experimental uncertainties 
in the future data scenario, we make the following as- 
sumptions: we assume that the systematic uncertainty 
on Aa^'p given in pQ improves by a factor ~ 2 and domi- 
nates the statistical uncertainty; furthermore, we assume 
several CP asymmetries to be me asured with the same 
precision as Aa^ P , see Table III The prospect for the 
uncertainty of the strong phase is taken from [67 . 

As in the AU = 1 model there is NP in the ss 
coupling, we assume that a^p(Z?° — > K + K~) and 
a^p(D + -> K S K + ) are enhanced, whereas a^p(_D° — !• 
7r + 7r~) and a,cp(D s — > Kstt + ) ~ 0. Furthermore, in 
the triplet model the observables a^ P (D + —> K$K + ) 1 
a^ P (D s -> K s tt+) and a^ p {D s -> K+ir°) are all deter- 
mined uniquely by A|, whereas does not contribute 
here. It is therefore difficult in the triplet model to 
account for these three CP asymmetries not being of 
the same order of magnitude. Consequently, we assume 
aQp(D s —> K + tt°) to stay large. All assumed central 
values are within the 2a interval of the current measure- 
ments listed in Table HT1 

Note that the expectation that the operator scus in the 
AU = 1 model contributes mainly to matrix elements 
with kaons is not a statement from SU (3) . It would 



amount to additional dynamical input, which we avoid 
in this work. Such assumptions could be implemented 
by assuming relations between matrix elements. 
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TABLE V: The coefficients Cd-ij of the physical S'C/(3)-breaking decomposition used in our analysis, see Section 11 B 
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